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Abstract 

Using first-principles calculations, the phonon frequencies at the T point and 
the dielectric tensor are determined and analysed for the cubic and rhombo- 
hedral phases of BaTiOa . The dipole-dipole interaction is then separated a la 
Cochran from the remaining short-range forces, in order to investigate their 
respective influence on lattice dynamics. This analysis highlights the deli- 
cate balance of forces leading to an unstable phonon in the cubic phase and 
demonstrates its extreme sensitivity to effective charge changes. Within our 
decomposition, the stabilization of the unstable mode in the rhombohedral 
phase or under isotropic pressure has a different origin. 

PACS numbers: 77.80.-e (Ferroelectricity and antiferroelectricity) , 
77.84.Dy (Niobates, titanates, tantalates, PZT ceramics, etc.), 63.10.+a (Lat- 
tice dynamics: general theory). 
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Barium Titanate (BaTiOs) is a typical ferroelectric material that undergoes three tem- 
perature phase transitions, from a paraelectric cubic phase, stable at high temperature, to 
ferroelectric phases of tetragonal, orthorhombic and rhombohedral symmetry. There have 
been considerable efforts to identify the origin of the transitions in this particular ABO3 
compound Among them, let us point out the seminal theory of Cochran In the 
framework of a shell model, he relates the ferroelectric transition to the softening of a trans- 
verse optic phonon at F. Within his model, the interatomic forces are decomposed into 
short-range forces and long-range Coulomb (dipole-dipole) interaction. The latter is evalu- 
ated by means of a Lorentz effective electric field, assuming a local spherical symmetry at 
each atomic site. Interestingly, the decomposition isolates the contribution of each kind of 
forces on the frequency of the transverse modes and identifies the structural instability with 
the cancellation of the two contributions. Although meaningful, Cochran's model is only 
qualitative. The numerical investigation is subject to many approximations. Moreover, it 
was shown by Slater |^ that the Lorentz field is far from spherical. 

In subsequent studies, it has been usually accepted that ferroelectricity in perovskites 
results from a delicate balance between short-range repulsions which favor the cubic phase 
and long-range electrostatic forces which favor the ferroelectric state. Although some cal- 
culations 0] illustrate this picture, they do not rely on a well defined separation of the 
interactions. 

In contrast, a separation of the interatomic forces was proposed recently by Gonze et 
al. p. Without postulating any atomic site symmetry, they introduce an analytic form 
for the dipole-dipole interaction at the microscopic level from Born effective charges and 
dielectric tensor. This formulation, evaluated thanks to first-principles data, generalizes 
Slater's calculation of the Lorentz field and can be used to refine Cochran's results 0. 

In our letter, we consider the cubic and rhombohedral phases as well as a compressed 
cubic structure. We first report dielectric tensor values, then compute the dynamical matri- 
ces and phonon frequencies at the F point and compare them with experiment. Separating 
a dipole-dipole interaction from the short-range remaining part of the dynamical matrix 
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following Gonze et al. p[, we quantify the balance of forces generating the unstable cubic 
phonon mode and investigate its sensitivity to effective charge changes. Within this decom- 
position, the hardening of the unstable phonon in the rhombohedral phase or under isotropic 
pressure has a different origin. 

Computational details are the same as in Ref. Calculations are performed within 

the Density Functional Theory (DFT) and the Local Density Approximation (LDA) using 
a conjugate- gradient plane-wave pseudopotential method. Responses to electric field and 
phonon-type perturbations are obtained within a variational approach to Density Functional 
Perturbation Theory. 

We consider cubic structures at the experimental and theoretically optimized volumes 
with lattice parameters = 4.00 and 3.94 A, as well as a compressed cubic cell with 
ao=3.67 A. For the rhombohedral phase, we worked at the experimental unit cell parameters, 
with relaxed atomic positions. Accurate values for the Born effective charges ^j! q,/? were 
already reported elsewhere [§. 

The computed dielectric tensor eoo is presented in Table | in comparison to its experimen- 
tal estimate in the cubic phase |0. No previous values were reported for the rhombohedral 
phase. It is well known that the DFT-LDA usually overestimates the experimental eoo- For 
the cubic geometry, the discrepancy is of the order of 25%. This error can be corrected in first 
approximation by the scissor oj»erator technique |jlO|. For all cases, we have used a scissor 



shift of 1.36 that adjusts the bandgap at the F point in the experimental cubic structure 
to the value of 3.2 eV For the cubic phase, the scissor corrected eoo (5.61) overestimates 
the experimental value (5.40) by less than 5%. For the rhombohedral structure, the values 
are globally smaller, especially along the ferroelectric axis. This goes hand-in-hand with the 
reduction also observed in this direction for Z* [^]. 

There are 12 optical phonons in BaTiOs. In the cubic phase, we have three triply 
degenerate modes of Fi^ symmetry and a triply degenerate mode of symmetry. Only 
the Flu niodes are infrared active, with LO-TO splitting, while the modes are silent 
modes that cannot be identified experimentally. Our values (Table H) at the optimized 



volume are in close agreement with the experiment [|TT| as the theoretical results of Zhong 



et al. |T2|. In particular, we reproduce the instability of the TOl mode corresponding to the 
vibration of Ti and Ba against the O atoms. The phonon frequencies appear very sensitive 
to the small volume change from the experimental to the theoretical cubic phase, contrary 
to Z* 1^ or 600- This is particularly true for the soft TOl mode, whose instability even 
disappears in our compressed cubic phase. 

Due to the long-range Coulomb interaction, the eigendisplacements of the TO modes 
(?7™) do not necessarily correspond to those of the LO modes ij]^^). The overlap matrix 
{rj^'~^\M\rj^^) reported in Table |T| establishes however that the mixing of modes is very 
weak in the cubic phase ( M is such that M = M^d^^/ with is the mass of atom k). 
In agreement with this observation, assuming that r^^*^ and rj^^ are identical, the fictitious 
LO frequencies predicted on the basis of the oscillator strengths (Eq. 10 of Ref. 0) are 
respectively of 701, 214 and 508 cm~^, in close agreement with the theoretical LO frequencies. 
Note the giant splitting of the TOl mode already mentioned by Zhong et al. [|T^. It arises 
from the large effective charges on Ti {Z^.- = +7.28) and O {Zq = —5.73, for a displacement 
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9.02. 



along the Ti-0 bond) generating a mode effective charge Z^q^ = — ^^toi|^toi) 

In the rhombohedral phase (Table each triply degenerate Fi^ (resp. F2u) mode from 
the cubic phase gives rise to a mode of Ai (resp. A2) symmetry with eigendisplacements 
along the ferroelectric direction, and a doubly degenerate mode of E symmetry. E and 
Al modes are infrared and Raman active. The only relevant comparative result we found 
is experimental |]13[ and localizes the phonon frequencies in three regions (100-300 cm~^. 



480-580 cm~^, and 680-750 cm~^) in qualitative agreement with our values. 

All the modes are stable in the rhombohedral structure. Due to the small distortions, 
the eigenvectors remain very similar to those of the cubic phase. Table fTT| compares Ai to 
corresponding Fi^ eigenvectors. Similar values are obtained for the E modes. They point 
out that Ai{T02) and E{T02) originate from the hardening of the soft mode. Even if both 
of these modes continue to couple strongly with the electric field, the smaller Z* make their 
mode effective charge smaller: 7.00 and 8.41 respectively. We predict a static dielectric 



constant equal to 33.09 along the ferroelectric direction and to 68.89 perpendicularly to it. 

The phonon frequencies uj and the associated eigendisplacements 77 are deduced from the 
dynamical matrix A through the following equation: Y.k'p -^ap^K, 1^') V^' ,P = M^u"^ ijf^^a- The 
a and (] indices denote the space direction while k and k' label the atom within the unit cell. 
The ansatz proposed by Gonze et al. can now be used to parametrize the dipole-dipole 
contribution to the interatomic force constant from the knowledge of and eoo, in the 



general case where these tensors are anisotropic |14 
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where Aq, = I]/3(£j»^)a/3 dp, d = Rj + Ti^i — 77, is the vector relating nuclei , and D = y A.d. 
The contribution of this dipole-dipole term to the dynamical matrix is evaluated using 
Ewald summation technique By this way, dipole-dipole (DD) and remaining short-range 
(SR) ||T5| parts of the dynamical matrix A can be isolated from each other d la Cochran 
{A = Add + ^sr) and their partial contribution to can be evaluated as follows: 

{r]\A\r]) = {r]\ADD\v) + {v\Asr\v\ . 

Add and Asr can then be modified independently in order to investigate their own influence 
on the instable mode. 

In Table |V| we report the values of oj'dd and for the TO modes of the cubic phase 
at the optimized volume. We observe that the instability of the FiuiTOl) mode originates 
from the compensation of two very large numbers, the DD interaction greatly destabilizing 
the crystal. Interestingly, this close compensation exists for the unstable mode only. 

In the cubic phase, the large values of Z^^ and (responsible of the strong Coulomb 
interaction) are mainly produced by a dynamic transfer of charge along the Ti-0 bond [Q. 
Postulating Asr to be fixed, we can fictitiously reduce this transfer of charge by decreasing 
simultaneously Z^^ and Zq^^, and monitor the Fiu{T01) mode frequency changes. Figure 1 
shows that a;^(T01) evolves quasi linearly with the transfer of charge. A change correspond- 
ing to a reduction of the order of 1% of Z^^ is enough to suppress the instability. Of course 



this situation is artificial: In a real material any modification of Z* would go hand-in-hand 
with a change of the SR forces. This result however highlights the very delicate compen- 
sation existing between DD and SR interactions. Interestingly, is also modified, due 
to the change of the eigenvector rj induced by the modification of A^d- This change is not 
crucial and a similar evolution of oj^ is observed if we keep the eigenvector of the original 
optimized structure. Note that all these conclusions are independent of the use of the scissor 
correction for eoo- From now, we report results without scissor correction. 

In the rhombohedral structure, there is no unstable mode although the eigenvectors 



remain close to those of the cubic phase (see Table It was found that Z* are smaller 
in this ferroelectric phase, suggesting a smaller DD interaction, but this could be partly 
compensated by a reduction of £00- For the Ai{T02) mode coming from the soft mode, 
(-286267 cm^^) add to a slightly larger SR counterpart (356373 cm^^). The values 
differ widely from those of the cubic phase: The SR forces give less stabilization but this is 
compensated by a larger reduction of the DD contribution. 

If we now fictively modify and replace Z* and of the ferroelectric structure by 
their value in the cubic phase, we modify the frequency of the Ai{T02) mode from 265 to 
266i cm~^: We obtain an instability even larger than in the cubic phase. From this point 
of view, the reduction of Z* in the rhombohedral phase appears as a crucial element to the 
stabilization of the Ai{T02) mode. Introducing Z* and e^o of the cubic phase, and 
are also strongly modified and becomes respectively -871017 and 800371 cm~^. The 
dramatic change of uj'g^^ results from a change of eigenvector pointing out the anisotropy 
of the SR forces (the overlap between the new and original eigenvector is equal to 0.86). 
If we had kept the eigenvector unchanged, we would still have observed an instability (74i 
cm^^) for the Ai{T02) mode although much smaller. This means that the inclusion of the 
effective charges of the cubic phase is already sufficient to destabilize the crystal but at the 
same time produces a change of eigenvector enlarging the instability. 

No more instability is present in the compressed cubic phase, although the global values 
of Z* do not differ significantly from those at the optimized volume f^. Moreover, the 



reduction of volume even increases the destabilizing effect of the DD interaction by 20%. 
However, strong modifications of the SR forces produce a mixing of modes so that no one can 
be identified with the unstable mode observed at the optimized volume. If we replace Asr 
by its value at the optimized volume we recover a very large instability (437i cm~^). The 
disappearance of the unstable mode under pressure seems therefore essentially connected to 
a modification of the SR forces in contrast to its stabilization in the rhombohedral phase 
related to the reduction of Z*. 
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FIGURES 

FIG. 1. Evolution of the Fi„(T01) mode frequency squared with respect to the dynamic trans- 
fer of charge along the Ti-0 bond (quantified here by the evolution of see text). Results are 
obtained with (o) or without (-•-) scissor shift for Cqo- SR and DD contributions to uJ^ (see text) 
are shown in the inset. 
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TABLES 

TABLE L Dielectric tensor of BaTiOs obtained within the local density approximation (LDA) 
or with an additional scissor correction (SCI). For the rhombohedral phase, the z axis points in 
the ferroelectric direction. 



Cubic phase [e^ 
ao = 3.67 A Oo = 3.94 A 



Rhombohedral phase 



LDA 
SCI 



6.60 
5.71 



6.66 
5.60 



= 4.00 A 
6.73 
5.61 



6.16 
5.26 



5.69 
4.91 



TABLE II. Phonon frequencies (cm ^) at the T point for cubic and rhombohedral BaTiOa. 



Cubic phase Rhombohedral phase 



Mode 


ao=3.67A 


Exp. Q 


ao=3.94A 


ao=4.00A 


Mode 




Mode 




Fiu{T01) 


214 


soft 


113i 


219i 


^i(roi) 


168 


£;(roi) 


161 


FiuiLOl) 


250 


180 


180 


159 


Ai{L01) 


180 


E{LOl) 


173 


Fiu{T02) 


296 


182 


184 


166 


Ai{T02) 


265 


E{T02) 


205 


Fiu{L02) 


513 


465 


460 


447 


Ai{L02) 


462 


E{L02) 


438 


Fiu{T03) 


737 


482 


481 


453 


Ai(r03) 


505 


E{T03) 


461 


Fiu{L03) 


1004 


710 


744 


696 


Ai(L03) 


702 


E{L03) 


725 


F2u 


308 


306 


288 


281 


A2 


274 


E 


293 
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TABLE III. Overlap matrix elements between the eigenvectors of the Fi„(rO) modes of the 
optimized cubic phase and those respectively of the associated Fi„(LO) mode and of the Ai{TO) 
mode of the rhombohedral phase. 



Fiu{L01) Fi„(L02) Fi„(L03) Ai{T01) A^{T02) Ai(T03) 



Fiu{TOl) 


0.17 


-0.36 


0.92 


0.13 


-0.97 


-0.19 


Fiu{T02) 


-0.99 


-0.07 


-0.16 


-0.99 


-0.13 


-0.01 


Fiu{T03) 


0.01 


-0.93 


0.37 


-0.02 


0.18 


-0.98 



TABLE IV. DD and SR contributions (see text) to the TO mode frequency squared (cm~^) 
for the cubic phase at the optimized volume. Values in brackets where obtained with the scis- 
sor-corrected value of Sno- 





Fiu{T01) 


Fiu{T02) 


Fiu{T02>) 


F2u 


^DD 


-625897 (-745610) 


7232 (8615) 


-130549 (-155518) 


109745 (130736) 


^SR 


613107 (732820) 


26538 (25155) 


361998 (386967) 


-26951 (-47942) 




-12790 


33770 


231449 


82794 
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